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Abstract 

The original Zee model can easily accommodate a bi-maximal mixing solution for the atmospheric 
and solar neutrino problems. From the most recent fit to the data we obtain a set of parameters 
for the Zee model. However, the recent claim on a positive signal observed in a neutrinoless double 
beta decay, which requires a nonzero mee, cannot be accounted for by the Zee model. Based on 
the observed neutrino mass-square differences and mixings, we derive a few general patterns for 
the neutrino mass matrix in the flavor basis, which can correctly describe the atmospheric, solar, 
and Oi^/3/3 data. Finally, we investigate a few possible extensions to the Zee model, which can give 
the required mass patterns. 
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I. INTRODUCTION 



More and more neutrino data are accumulated since a few years ago. They all contribute 
to the understanding of neutrino masses and mixings. They are very important clues to the 
underlying theory for neutrino masses, and perhaps to other fermions as well. It is therefore 
very important to pin down the necessary patterns for neutrino mass matrix that can explain 
all the observations. 

Deficits in both atmospheric and solar neutrinos have been known for a long time. A 
possible and favorable explanation is neutrino oscillation, which is made possible because 
neutrinos have different masses and the flavor basis is not the same as the mass basis. The 
atmospheric neutrino data from SuperK [Q and the preliminary data from K2K [Q showed 
a maximal mixing between the and Uj- with a mass-square difference and a mixing 
angle: 

Am^tm ~ 3.0 X 10-3 eV^ , sin^ 2e^t^ = 0.99 . (1) 

On the other hand, all the combined solar neutrino data favor the oscillation of the z/g into 
a mixture of and u^-. Especially, the recent SNO |Q data provided a convincing evidence 
that there is a non-z/e component in the neutrino flux coming from the sun. The most recent 
fit favors, at least statistically, the large mixing angle (LMA) solution, namely, [^] 

Amj^i ^ 4.5 X 10"^ eV^ ; tan^ ^^^i = 0.41 , (2) 

which gives sin^ 2^soi ~ 0.8. The above two mass-square differences can be accommodated 
by 3 active neutrinos. The much debated result from LSND if confirmed, requires an 
additional neutrino, which is sterile. 

Another recent result from a neutrinoless double beta decay (Oz//3/3) experiment added 
another clue or constraint to the neutrino mass matrix. The positive signal in the Oi>f3(3 decay 
implies a nonzero ruee entry in the neutrino mass matrix arranged in the flavor basis: 

^ '^ee ^e/^ ^er 
^efi 

\ f^er ^^/iT fTT'TT 

The 95% allowed range of rriee is 0.05 — 0.86 eV with a best value of 0.4 eV. It also implied 
that the electron neutrino is majorana in nature. There is, however, an argument against 
this claim: see Ref. [|, H. 



Such a large value for rriee ^ yAml^^^^^^, if true, gives a nontrivial modification to the 
neutrino mass patterns and mixings. A lot of possible mass textures that were proposed 
to explain the atmospheric and solar neutrino deficits become incompatible with the new 



OuPP data. A well-known example is the Zee model |T0[ that can generate a bi-maximal 
mixing between the Un and Uj. and between the and the mixture of z/^ and [O, |12|, |T3f ^ 



however, the Zee model guarantees the diagonal mass matrix elements to be zero. Therefore, 



^ For a number of analyses on the Zee model after the SNO data but before the 0i//3/3 data, please see Refs. 

H H, @ 
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it cannot explain the nonzero rriee implied by the 0^/3(3 signal. It would then be necessary 
to consider additional new physics to the Zee model in order to accommodate the uiee data. 
In this work, we look into the general mass patterns that can accommodate all the data, 
and from the patterns we introduce additional new physics to the Zee model. 

The organization is as follows. In Sec. II, we describe the Zee model and how it could 
explain the atmospheric and solar neutrino data. Using the most recent fit to the data we 
obtain a set of values for the parameters of the Zee model. In Sec. Ill, we derive the general 
mass patterns that can explain all the data. In Sec. IV, we point out a few simple extensions 
to the Zee model that can generate the necessary patterns. We conclude in Sec. V. 



II. ZEE MODEL 



An economical way to generate small neutrino masses with a phenomenologically favor- 

% O, which generates masses via one- loop 



able texture is given by the Zee model |TT 



diagrams. The model consists of a charged gauge singlet scalar h^ee, the Zee scalar, which 
couples to lepton doublets ipL via the interaction 



fab (V'^lC^: 



/i7,,, +h.c. 



"'Zee 



(3) 



where i, j are the SU(2) indexes, a, h are the generation indexes, C is the charge-conjugation 
matrix, and fab are Yukawa couplings antisymmetric in a and h. Another ingredient of the 
Zee model is an extra Higgs doublet (in addition to the one that gives masses to charged 
leptons) that develops a vacuum expectation value (VEV) and thus provides mass mixing 
between the charged Higgs boson and the Zee scalar boson. Let us parametrize this mixing 
term as 

K Kee + h.C. , (4) 

where is of dimension [mass]. This coupling, together with the /ab's, enforces lepton- 
number violation. Through one-loop diagrams the Zee model can generate off-diagonal 
majorana mass terms, mab{ci, b 



e, /U, r), given by 0, [T|, [19 



1 



mab 



167r2 y/2Mwcos(3 



gml 



OHM' Ml- 



) 



(5) 



where the function h{x,y) = \og{y/x)/{x — y) and the charged Higgs h'^ couples to the 
charged lepton / with a coupling grrii/ {\/2Mw cos /?). It is a generic feature of the Zee model 
that the diagonal terms rriaa = 0. Thus, we obtain a neutrino mass matrix of the following 
form 



/ 



m. 







m, 
m 



\ me 



m 



Recent analyses [TT|, 
of the following texture 



14, 15, 16, 17,^ 



211 showed that the Zee mass matrix 



/ 



m 





e 



e 
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where e is small compared with me^^ and mgr, is able to provide a compatible mass pattern 

is required to 



that explains the atmospheric and solar neutrino data. Moreover, 
give the maximal mixing solution for the atmospheric neutrinos. 

Keeping the leading power in e we obtain the mass eigenvalues mi, 1712, and the diag- 
onalization matrix U such that U^M^U = diag{mi,m2,ms} (we work in the basis where 
the charged- lepton mass matrix is diagonal): 



7712 

U 



m 



m 



m 



er 
irier 



-2e 



( *(1 



2m3 



rrier 



+ e 
+ e 



mer{2ml^-ml^) me 



^ ) 

er / 



V2n 



+ e 
+ e 



2m3 



2\/2m 



2 ) 

er / 



J 



The above results reduce to the results of Ref. [|13] in the limit 



where m = 
e 0. 

The mass eigenvalues mi, m2, and U are arranged in the inverted mass hierarchy such 



that mf 



^ 1712 "^3- 



The mass-square difference between mf and rn^ explains the solar 



neutrino while that between rn^ and m| explains the atmospheric neutrino. The mixing 
angles can also be obtained from the matrix U such that tan6'soi = tan ^12 = Ue2/Uei and 
tan6'atm = tan 6^23 = U^^/Ur^. According to the most recent fit to the atmospheric neutrino 
data in Eq. (P, the sin^ 26'atm is practically 1 and we obtain 



m„ 



±1 and 



m 



m 



e/i I 



m„ 



0.039 eV . 



(6) 



The condition |m 



en I 



mp 



also implies t/es = 0, which is consistent with the CHOOZ limit 



2^ . The parameter e is then given by 



m, 



e/i 



4.1 X 10"^ eV for LMA 

4.2 X 10-9 eV for VAC 



(7) 



We note that the mixing angle tan6'i2 = 1 + e/(2^/^me^) essentially equals 1, which is 
obviously incompatible with the SMA solution, but fairly consistent with the LMA and 
VAC solutions Since the LMA solution is the more favorable one, at least statistically, 
and moreover the required e for the VAC solution is a much more fine-tuned one, we shall 
concentrate our discussion on the LMA solution in the next section. 



We can use Eq. (||) to obtain the values on fat from the fitted values of 



m, 



The ratio mg^ : mg,- : m^^ 



1:1: AmgQi/(\/2Am^tjjj), which implies 



m 



fell ■ fer • fi 



m. 



sol 



tJ.T 



ml ml A/2Am^ 
1 : 3.5 X 10~^ 



atm 



3.7 X 10~^ (3.8 X 10 



-10\ 



(8) 
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where the numerical value in the parenthesis is for the VAC solution. Taking 
Mf^-,M,- ,m3 ~ 100 GeV, we obtain /e^ ~ 10"^ - 10-^ fer ~ 10"^ and /^^ ~ 10^^ 

d Zee 

for the LMA solution. 

These fabS are also subject to other constraints coming muon and tau decays ||23[. The 
decay ^ — > euu puts a bound on /g^: |/e^p/M^ < 7 x 10~^Gf, where M is a common 
mass scale for M^-, M^- and Gp is the Fermi constant. Taking M ~ 100 GeV, we obtain 

d Zee 

\fefi\ < 9 X 10^'^, which is consistent with the /e^ obtained above to explain the neutrino 
mass. The branching ratio B{fi — > 67) puts a bound on |/^T-/eT-|/M^ < 2.8 x 10~'^Gf, which 
implies |/^r/er| < 3 X 10~^ for M ~ 100 GeV. This is also consistent with the f^T,fer 
obtained above. 



III. PATTERNS PREDICTED BY Ouf3f3 DATA 



If the 0z//3/3 signal is true, it demands a very large value for niee relative to the mass- 
square differences required for the atmospheric and solar neutrino data. It has been shown 



2J, ^ |2^, ^ that the most favorable mass pattern is a degenerate one with small mass 
differences among the mass eigenstates. ^ 

To obtain the mass pattern it is easy to work bottom-up. We work in the basis where 
the charged-lepton mass matrix is diagonal such that the mixing information is entirely 
contained in the neutrino mass matrix. From the data we can write down specific forms 
of the mixing matrix U and the diagonal mass matrix in the mass basis, then we can 
obtain the mass matrix Mi, in the flavor basis as 



M,, = U Md 



where the flavor basis and mass basis are related by 



(9) 



Ml 



f/e2 








S) 




f/r2 




\ V3 



(10) 



Assuming no CP violation a general form of U is given by 

U = U23 Ui3 Uu , 



(11) 



where Uij is the rotation matrix about the i and j'th mass eigenstates. By considering the 
various data sets the general form of U can be reduced to a simple form. Since we know that 
the atmospheric neutrino requires a maximal mixing between the and u^, U23 is given by 



23 



[10 \ 
J L 

^ V2 V2 

^ ^ 



^ To explain the Qvj3l3 data alone it is not necessary to invoke neutrino oscillations. For example, i?-parity 
violation M and other new physics such as leptoquarks [M can explain the QvPP data. 
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Since the angle ^13 is constrained by CHOOZ |22[ to be small, we simply set Ui^ = I. For 
the rotation between the (1, 2) states there are a few possible solutions to the solar neutrino, 
and so we used a generic Uu as 



U- 



12 



/ c s 0\ 

-s c 

V 1/ 



where c = cos0i2, s = sin 6*12. Therefore, U is given by 



U 



I 



\ 



s/2 

s c 
' V2 V2 



\ 

1 

73 ) 



(12) 



For the mass matrix in the mass basis there are two cases: (i) normal and (ii) inverted 
mass hierarchies, which we shall consider in turns. In the normal hierarchy, the solar os- 
cillation is between the two lighter states, while the atmospheric oscillation is between the 
heaviest and the lighter states. Recall the convention that the solar oscillation is between 
the "1" and "2" states, we put the diagonal mass matrix as 



j^^normal 



mo 

(mo + 6) e'^ 

(mo + m) e'^' 



(13) 



where 5 



Amii/(2mo) 



Am^tjn/ (2mo), and for cj) and 



We shall later show the results for other values of (f) and 
we obtain the neutrino mass matrix in the flavor basis: 



i)' we take as for the moment. 
Using Eqs. (p 



M, 



normal 



\ 



With input values of mo 
0.41/1.41, we obtain 



f2^ 

0.4 eV, 6 



2 ' 2 
+ f + J 



T2|), and (0), 



(14) 



5.6 X 10-5 eV, m = 3.75 x 10^^ ^y, c = 1/^141, s 



Am 



32 



Am 



21 



= 3 X 10"^ eV^ 
4.5 X 10-^ eV^ 



sin^ 2^23 
tan^ 9i2 - 



-- 1 
0.41 



which are exactly the best fits shown in Eqs. (0) and (0). This is obvious because the input 
values for the parameters that we chose are simply based on the data. In other words, with 
the data we can fit the parameters of the model. 

For the case of inverted mass hierarchy the solar neutrino oscillation is between the two 
heavier states while the atmospheric neutrino oscillation is between the lightest and the 
heavier states. Recall the convention again that the solar oscillation is always between the 
"1" and "2" 



states, we used a M|^™=^*^'^: 



M 



inverted 
D 



/ mo + m + 6 
mo + m 

V mo 



(15) 
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in which we have chosen the phase angles to be zero. Using Eqs. 
obtain 



dll), and (|T5D we 



M, 



inverted 



/ 1710 + 171 + 6 --^ 5 



— S 



mo + f + 



-— 8 

'^5 ^ + 45 



With the same input values for mo, 6, m, c, s we obtain 



'''32 



Am"- ~ 



Am^i 



-3 X 10"^ eV" 



sin^ 2^23 



Am^i ^ -4.5 X 10"^ eV^ , 



tan^ 012 = 0.41 

In general, we can use the data to fit the values of mo, m, 6, 612: 

mo = mee ; m ~ Aml^^/ (2mo) ; 
6 ~ Am2^i/(2mo) , 9^ = ^soi • 



(16) 



(17) 



We have already used a particular form of U such that a maximal mixing is always between 
the "2" and "3" states. 

Next we are going to show the results when we allow nonzero phases for (p and </)'. For 
illustrations we use (ii) = vr, 0' = 0, (iii) = 0, 0' = vr, and (iv) = vr, 0' = vr. The case 
(i) = 0' = has already been shown above. The neutrino mass matrices in the flavor basis 
for normal and inverted mass hierarchies are given, respectively, by 
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m 



(iv) 



= vr, 0' = : 



cos 26'i2mo 

sin 2612 



■5 3" 



sin 29i2 ^ cs s: 

mo - ^ d 



cs s: sin20i2^ cs s: \ 

72^ 72^"^0~72^ \ 

V2 -| COS 2^12^0 + ^2o±Hi_|5 _i COS 2^i2mo-^2o±i£_ ^5 

-^ii^mo -^S -i COS 2^i2mo - - cos 2^12^0 + - 
cos2^i2(mo + m) + 5c2 --^n^^jrio + m) - ^6 
-^ii^(mo + m) - ^(5 -i cos 2^12 (mo + m 

sin 2612 



V2 

= 0,0' = 
mo + (5 



V " ■ ' V2 -2 

(mo + m)-^5 -i cos 26112 (mo + 



/ 

'-^{mo + m)-^S 
+ 4''+T'^ -icos2^i2(mo + m)- ™' 



m) - ^ + 



2 ^ 2 — " V--U ' ■■ V 2 ' 2 " 

45 COS 20i2(mo + m) + ^ + ^6 



+ ^6 



IT : 



-f + f (5 mo + 



mo + f + f(5 -f + f^ 
mo + m + 5c^ 



— — 



/ 

= TT 



m _|_ s f 
2 2 " 

mo + f + 



mo + f + %5 



^ + ^6 



J 



' = TC,(b = 



I cos 26'i2mo — 5s 

1 siii2ei2 cs s: 1 

V 72^"^0 ^72^ ~ 2 



sin 2012 ^ cs s: 

^ii^mo - ^ 5 -i cos 2^i2mo - - y ^ -| cos 2^i2mo + 

sin 201 2 ^5 _icos2^i2mo + - ^5 -icos2^i2mo - 



_ sin 2012 



a 2012 _ cs r \ 
72^"^0 72 ^ \ 
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-^ii^(mo + m) - ^5 -lcos29u{mo + m) - ^ + ^6 -i cos 2^12(^0 + m) + ^ + ^5 
^-5Hi|ii(mo + m) - ^5 -icos2ei2(mo + m) + ^ + 45 cos2^^i2(mo + m) - + ^5 y 

The present data cannot distinguish the above patterns because they all give correct 
Am^^jjj, sin^26'atm, ^"^soi; tan^ 6'soi, and rriee by adjusting the parameters mQ,m,6,9i2- The 
preference of one pattern over the others is a matter of taste. We suggest two fine-tuning 
criteria: (1) the smaller the ratios m/mo and 6/mQ the worse the fine-tuning is, and (2) 
the stability of the maximal mixing provided by the mass matrix. The first criterion prefers 
cases (i) and (iii) because cases (ii) and (iv) require a larger mo in order to explain the Oz//3/3 
data. Second, the mass matrices in (i) are not stable, which is explained as follows. Let 

us look at the (2,3) block of the matrix, which is of the form ^ ^, where e <^ 1. The 

maximal mixing between the "2" and "3" states is very unstable. Once m^^ and differ 
very slightly, the maximal mixing is destroyed. We have verified that if m^^ or rrirT increases 
by 1% the resulting sin^ 26'atm ~ 0.5, already out of the allowed range. It is barely consistent 
with the data range if m^^ or rrirT is changed by less than 0.5%. On the other hand, in 

case (iii) the (2, 3) block is of the form ^ ^ This is a rather robust form such that even 

when m^^ or m^-T- is changed appreciably, the maximal mixing remains. We have checked 
even if m^^ or 171^-^ increases by a factor of two the maximal mixing solution remains for 
the atmospheric neutrino. This is easy to understand. As long as m^^ ^ ?7i/j/i,^rT we have 
the maximal mixing in the (2, 3) block. So we conclude this section by saying that the mass 
matrices in case (iii) are the more favorable ones in view of stability and fine tuning. 

Since in order to explain the Oi>f]f] data a nonzero rriee is needed, which immediately rules 
out the original Zee model. In the next section, we look at a few extensions that can explain 
the Oz//3/5 data while maintaining the solutions to atmospheric and solar neutrino problems. 
We end this section by saying that we have checked the patterns given in Ref. We 



found that we can reproduce their Aml^^, Amg^j, and sin^ 26'atm, but not sin^ 26*801- We 
believe they calculated the solar mixing angle using the wrong pair of states, instead of 
using the pair whose mass-square difference is Amg^j. 



IV. EXTENSIONS TO THE ZEE MODEL 

Since the contributions to neutrino mass matrix in the Zee model arise from 1-loop 
diagrams, it is natural to consider the 2-loop contributions as well. It was pointed out 



that some 2-loop generalizations of the Zee model may be able to generate nonzero diagonal 
matrix elements. However, it is very difficult to accommodate the LMA solution unless by 



fine tuning of the model parameters [^]. In the following, we consider a couple of extensions 
that can hopefully accommodate the data in a less unnaturally way, though fine tuning may 
still be necessary. 

A straight-forward extension is to add a heavy right-handed neutrino Ni for each gen- 
eration with an interaction: C = —yiLiHuNi — ^NiNi + h.c, where Mj is a mass scale 
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much larger than the weak scale and i/i is a Yukawa coupling. After the electroweak sym- 
metry breaking, the interaction gives a Dirac mass term to the corresponding SM neutrino. 
Through the see-saw mechanism a small neutrino mass is then generated. Integrating out 
the heavy right-handed neutrinos, the effective interaction term that gives rise to neutrino 
mass is given by 



2 {LiHu){LiHu) 



C = y T' l^;^'^^""' + h.c. , (18) 



where we shall simply take a common scale for all Mj = M. In this way diagonal majorana 
mass terms can be generated easily as 

m.. = -y^^ , (19) 

for i = e, fi, T. The appropriate values for the diagonal matrix elements displayed in the 
last section can be obtained by adjusting the ?/j's, though in some sense it may require a 
fine-tuning. The off-diagonal matrix elements can be adjusted by varying the parameters of 
the Zee model. 

Another interesting extension to the Zee model is in fact a supersymmetrization of the 
Zee model within the minimal supersymmetric model (MSSM) without i?-parity conserva- 
tion [|l^. The right-handed slepton works as the charged gauge singlet (Zee singlet) while 



the MSSM already has two Higgs doublets, which provide the mixing between the and 
Zee singlet through the VEV of the Here the i?-parity violating terms fiiLiH^ provide 
mixings among neutrinos and higgsinos. Through a see-saw type mechanism tree-level con- 
tributions to neutrino masses arise and the rriij are roughly proportional to Hi/ij/ fi, where /i 
is the Higgs parameter in the ixHuHd term. Thus, the diagonal matrix elements are nonzero. 
There are also one-loop type contributions to the matrix elements in both diagonal and off- 
diagonal elements. It is then possible to adjust the soft and i?-parity violating parameters in 
order to obtain the mass matrix in (iii) of the last section. Again, to some extent fine-tuning 
of the parameters is needed to ensure the form of the mass matrix. 



V. CONCLUSIONS 

We have examined the Zee model in light of the most recent fits to the solar and atmo- 
spheric neutrino data. The Zee model is a natural candidate that can provide a bi-maximal 
mixing solution to the atmospheric and solar neutrino problems. The resulting atmospheric 
mixing angle ^23 is right at the best value of the experimental data while the solar mixing an- 
gle 9i2 is fairly consistent with the allowed range of the LMA and VAC solutions. By fitting 
to the data we found that the lepton-number-violating coupling /e^ is of order 10^^ — 10~^ 
and the ratio fe^i : /er : /^r = 1 : 3.5 x 10"^ : 3.7 x 10"^ for the LMA solution. The size 
of these couplings is consistent with other constraints coming from muon and tau decays. 

The recent neutrinoless double beta decay 0, ||, ||, if confirmed, rules out the original 
Zee model, because the positive signal requires a nonzero niee matrix element. Moreover, a 
large rriee value necessarily gives large mi^2,3 for the mass eigenstates. Such neutrino masses 
form a significant fraction of the warm dark matter |3^. We have derived a few general 
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patterns of neutrino mass matrix from the observed data of neutrino oscillations (here we 
assumed the LMA solution for the solar neutrino.) Due to fine-tuning and stability reasons 
the forms in (iii) of Sec. Ill are preferred. Nevertheless, the interpretation of the Oz//3/5 
signal requires a consensus among various experimental groups before generally accepted by 
the community. 

Finally, we offered a couple of extensions to the Zee model that can give nonzero diagonal 
matrix elements. However, in order to obtain the appropriate mass patterns fine-tuning of 
the parameters involved is still necessary. 

We are grateful to Otto Kong for useful discussions. This research was supported in part 
by the National Center for Theoretical Science under a grant from the National Science 
Council of Taiwan R.O.C. 
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